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Abstract 

One of the interesting properties of the NMSSM is that it can accommodate a light pscu- 
Q\ ' doscalar of order 10 GeV. However, such scenarios are challenged by several experimental con- 

straints, especially those related to the fermionic decays of the pseudoscalar. In this Letter, 
we extend the NMSSM field content by two gauge singlets, with lepton numbers +1 and —1. 
. Lhis serves the twin purpose of generating neutrino masses via the inverse seesaw mechanism 

Oh' an d keeping the option of a very light pseudoscalar experimentally viable by opening dominant 

(JL) i invisible decay channels of the pseudoscalar which help it evade the existing bounds. 

<d : 

^ . PACS Nos: 12.60.Jv, 14.80.Cp, 14. 60. St 
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We consider the Next to Minimal Supersymmetric Standard Model (NMSSM) which provides a 
natural solution to the so-called \i problem through the introduction of a new gauge singlet super- 
field S in the superpotential. In the NMSSM, the \x parameter is linked to the vacuum expectation 
\ value (VEV) of the scalar component of S whose size is of the order of the supersymmetry breaking 
scale pQ. Since it permits a scale invariant superpotential, the NMSSM is the simplest supersym- 
metric generalization of the Standard Model (SM) in which the supersymmetry breaking scale is 
the only mass scale in the Lagrangian. Moreover, by providing an additional tree level contribution 
to the quartic term of the scalar potential it can ameliorate the 'little hierarchy problem' of the 
MSSM, which is related to the requirement of a large (^> Mz) soft supersymmetry breaking mass 
that can push the mass of the lightest neutral Higgs beyond the LEP-2 limit of 114 GeV [2]. Fur- 
thermore, the NMSSM can admit a very light CP-odd Higgs boson (m Al ~ 1 - 10 GeV) [31 E]. The 
main objective of this Letter is to improve the experimental viability of such a light pseudoscalar 
by providing it with a dominant invisible decay mode in a minimal extension of the NMSSM which 
contains a source of lepton number violation that also yields an acceptable neutrino mass. 

In the NMSSM, the lightest CP-odd physical scalar A\ can be decomposed as 

Ai = cos 6 A A M SSM + sin 9 A A S , (1) 

where ^4mssm is the MSSM part of the CP-odd scalar, which arises solely from the NMSSM Higgs 
doublets, and As is the part that arises from the new singlet superfield S. It is the singlet admixture, 
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i.e. the sin^ projection, that allows the NMSSM pseudoscalar to be much lighter than what it 
could have been in the MSSM. On the other hand, if Ai is very light then its detection crucially 
depends on its couplings to quarks and leptons, which depend on cos 6^. These couplings can be 
extracted from the following part of the Lagrangian [5] : 



C A 1 ff= X u(d)j^fl5fA 1 , (2) 



where g is the SU(2) gauge coupling, Xd(X u ) = cos 8 a tan j3 (cos 8 a cot j3) for down-type (up-type) 
fermions, tan/3 = v u /vd with v u and denoting the up- and down-type Higgs VEVs. 

A light pseudoscalar is phenomenologically interesting mainly for two reasons: 

(i) The direct search limit of 114 GeV on the mass of the SM-like Higgs (a slightly smaller limit of 
93 GeV for the MSSM neutral Higgs) is obtained from its non-observation at the highest energy 
run at LEP-2, where the Higgs was expected to be produced from gauge interactions (full strength 
ZZh coupling) and decay dominantly (~ 75%) into bb final states. In the NMSSM, there could be 
two important changes: 

• The lightest CP-even Higgs (h) may have a large singlet component, which leads to a signifi- 
cant dilution of the ZZh coupling. The Higgs production cross section will then be reduced, 
and hence the lower limit on rrih will be relaxed. 

• h may dominantly decay into a pair of A\, with each A\ decaying into //, where / is b, 
t, c, or /i, depending on the kinematic thresholds. Therefore, the existing LEP-2 Higgs 
search strategy from 2b final states would fail as one should look for 4/ final states. In fact, 
a reanalysis of the LEP-2 data by the LEP Collaboration has already put constraints on 
o°m{<^ <^->Zh) x Br(/i ~~ -^l-^i) x Br(Ai — > ff) 2 , where / = b or r depending on kinematics 
thresholds [H [7]- Similarly, for m J 4 1 < 2m T , upper limits have been placed on a(pp — >• 
hX) x Br(/i — > AiAi) x Br(yli — s> fi + /i~) 2 by the DO collaboration at Fermilab Tevatron [8]. 
The upshot is that the lower limit on rrih is slightly weakened, and a smaller value of rrih (e.g. 
105 GeV) is in fact preferred by electroweak precision tests [H] ■ 

(ii) If the lightest supersymmetric particle (LSP) happens to be very light (a few GeV), then a 
light A\ offers the possibility of s-channel LSP pair-annihilation into an on-shell A\. This resonance 
channel has a special significance when one attempts to account for the observed dark matter relic 
abundance. It has recently been shown that a light LSP of mass ~ 10 GeV can have interesting 
consequences in the context of the recent DAMA/CoGeNT results |10j . 

Let us now briefly discuss the existing bounds on the mass of the light pseudoscalar. The 
constraints on X^, defined in Eq. ([2]), for approximately in the range of 1 to 10 GeV have been 
summarized in [HE]. Measurements of AM^ S , Br(£? — > X^), Br(l? + — > t + v t ), and particularly, 
Br(i3 s —> fi + fi~) severely constrain [H]. The rates of these processes primarily depend on 
the choice of tan f3 and the soft supersymmetry breaking trilinear term At, and the constraints are 
in general weaker when these parameters are small. Values of itla 1 between 1 GeV and m;, are 
generally disfavored from £>-meson data [1] . Constraints on mA± also arise from radiative T decays 
P21 [13l [H] , namely, T(nS) — > 7A1, with A\ — > /i + /i~(r + T _ ) (further investigated and reviewed 
in [5]). Severe constraints also arise as a consequence of rjb — A\ mixing [151 [T6| I17j. The different 
mA 1 windows which are sensitive to different processes are listed in Table [TJ The table also shows 
the ranges where the LEP (ALEPH [7] and OPAL |18j ) constraints are applicable. The origin of 
all these constraints can be traced to the visible decay modes of A\. 
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Processes 


m^j < 2m T 


[2m T ,9.2 GcV] 


[9.2 GeV,A/ T(ls) ] 


[M r(ls) ,2m B ] 


T jAi -> 7 + (u + u~,gg, ss) 


/ 


X 


X 


X 


T — > jAi — > 7t + t~ 


X 


/ 


X 


X 


yli-?yb mixing 


X 


X 


/ 


X 


e+e- -» Z + 4r (ALEPH) 


X 


/ 


/ 


/ 


e+e- -> 66t+t- (OPAL) 


X 


X 


/ 


/ 



Table 1: Different processes constraining different m,A 1 windows. The V" symbol in a given entry attests 
the existence of important or meaningful constraints from a given process, while the "x " symbol implies 
otherwise. 



However, the situation may dramatically change if A\ has dominant invisible decay modes. Its 
decay into a pair of stable neutralinos (if kinematically possible) is one such example. The BABAR 
Collaboration [HJ] at the PEP-II 5-factory has, however, searched for radiative T-decays where 
a large missing mass is accompanied by a monochromatic photon, and from its non-observation 
has set a (preliminary) 90% C.L. upper limit on Br(T(3S') — > 7^4i) x Br(Ai — > invisible) at 
(0.7 - 31) x 10~ 6 for m Al in the range of 3 to 7.8 GeV. 

In this Letter we further explore the possibility of invisible decay channels that would allow a light 
Ai escape detection even outside the range of 3 to 7.8 GeV. We show that if we extend the NMSSM 
by two additional gauge singlets with non-vanishing lepton numbers, they would not only provide 
a substantial invisible decay channel of A\ but, as a bonus, would also generate small neutrino 
masses through lepton number violating (AL = 2) interactions. The visible decay branching ratios 
of Ai would then be reduced. As a result, the constraints on would be weakened. A light A\ 
can then be comfortably accommodated. 

In the framework of the NMSSM, neutrino masses and mixings can be generated via different 
mechanisms, either under the assumption of i?-parity conservation or i?-parity violation (RpV). 
In the latter case, light neutrino masses/mixing data can be successfully accommodated via the 
inclusion of explicit trilinear and/or bilinear RpV terms [20], or through spontaneous violation 
of lepton number in the presence of right-handed neutrino superfields (in addition to the NMSSM 
singlet S) [21] . Alternatively, one can employ a standard seesaw mechanism in a -R-parity conserving 
setup by adding three gauge singlet neutrino superfields Ni to the NMSSM particle content [22]. 
In this case, the light neutrino masses originate from the tiny Yukawa couplings and (dynamically 
generated) TeV scale Major ana mass terms. 

Here, we follow none of the above paths. We rather implement the "inverse seesaw" mechanism 
[23] by adding two gauge singlet superfields with opposite lepton numbers (+1 and —1). During 
this implementation we assume that i?-parity is conserved, i.e. we do not admit any AL = 1 term 
in the superpotential. To appreciate the advantages of the inverse seesaw over the standard one we 
look at the difference, from the point of view of effective operators, between the properties of the 
dimension-5 (d-5) Weinberg operator which is lepton number violating and the dimension-6 (d-6) 
operator which is lepton number conserving. In the d-5 case, the light neutrino mass is given by the 
seesaw formula m u ~ m^/M, where mr) = fv is a Dirac mass, with v as the electroweak VEV and 
/ as a generic Yukawa coupling. The source of lepton number violation is the Majorana mass M. If 
we demand the Yukawa coupling / to be order one, then for m v to be ~ 1 eV, M has to be close to 
the gauge coupling unification scale. On the other hand, the d-6 operator goes as 1/M 2 , but since 
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this operator conserves lepton number, its coefficient is not related to that of the d-5 operator |24| . 
Therefore, for the d-6 case, one needs a separate source of lepton number violation to generate light 
neutrino mass. The generic form of this mass is m v ~ (m^/M 2 )^, where fi is a lepton number 
violating (AL = 2) mass parameter. In this case, one can comfortably keep the fundamental scale 
M close to TeV, i.e. within the LHC reach, and yet choose / to be order unity which can trigger 
large lepton flavor violation. Here, the lightness of the neutrino mass (eV, or even lighter) is due to 
the smallness of fi. Having this small dimensionful term in the Lagrangian is technically natural in 
the sense of 't Hooft [25], as in the limit of vanishing fi one recovers the lepton number symmetry. 
The structure of light neutrino mass obtained via inverse seesaw conforms to the principle of 'low 
fundamental scale, large Yukawa coupling and light neutrino mass' [26j . 

In our model, the superpotential is given by 

W = Wnmssm + W , where, (3) 

Wnmssm = ffjHdQiDj + fljHuQiUj + fljR d LiEj + \ H SH d H u + -S 3 , (4) 
W = fijH u LiNj + (XN)iSNiXi + nxiXiXi . (5) 

In the above expressions, H d and H u are the down- and up-type Higgs superfields; Qi and Lj 
denote the SU(2) doublet quark and lepton superfields; Ui (D{) and Ei are the SU(2) singlet up 
(down)-type quark superfields and the charged lepton superfields, respectively. We have denoted 
the Yukawa couplings by / with appropriate flavor (u, d, e, v) and generation indices = 1, 2, 3). 
S is the singlet superfield already present in the minimal NMSSM. Besides, we have added two more 
gauge singlets N and X, for each generation, which carry lepton numbers L = — 1 and L = +1, 
respectively. In our formulation, even though L is not a good quantum number because of the 
presence of a non- vanishing fix > ( — 1) is still a good symmetry. We have written the NX and XX 
terms in a generation diagonal basis without any loss of generality. Once the scalar component of 
S acquires a VEV (vs), not only the conventional fi-term is generated with fi = Xh^s, a lepton 
number conserving mass term Mx^n^x is generated as well, with Mx = Xxvs- One more lepton 
number conserving mass term mo^u^N emerges with mr> = f v v u . 

The crucial term relevant for inverse seesaw is the AL = 2 term involving fix, which is the only 
mass dimensional term in the superpotential. We assume that the Z3 symmetry of the superpo- 
tential is absent only in this term. We treat fix as an extremely tiny effective mass parameter 
generated by some unknown dynamics. Its smallness would eventually decide the lightness of the 
light neutrino. We make a few observations at this stage: 

(i) The superpotential treats the two singlets N and X differently in the sense that it yields 
a [ix^ x^x Majorana mass term (AL = 2) but does not lead to a similar fix^N^N term. This 
discrimination requires further qualification. A generic superpotential with (— 1) L parity should 
have included the latter term. Both fix and fix can be naturally small, as their absence enhances 
the symmetry of the Lagrangian. But the important thing to note is that the magnitude of fix 
(and not that of fix) controls the size of the light neutrino mass \27\ 128] , In view of this, for the 
sake of simplicity, we have assumed fix = 0. 

(ii) Two questions naturally arise here. First, although we have put fix by hand and claimed 
that it is tiny, is it possible to dynamically generate a small fix starting from a superpotential which 
does not a priori contain any mass dimensional term? Second, is it possible to realize fix *C fix in 
a sensible model? 



4 



To address the first question, we admit that with the particle content of our model it is not 
possible to provide a natural solution for a small fix- One possibility could have been to start with 
a trilinear XSXX term in the superpotential of Eq. (|5|), which would lead to fix = ^ v s- Since 
v$ ~ v, the requirement to produce the correct light neutrino mass would then compel us to take 
A ~ 10 Such a small dimensionless trilinear coupling would be against the spirit of inverse 
seesaw mechanism as illustrated earlier. However, by extending the particle content of the model it 
is possible to dynamically generate a small fix- 111 this context, we recall that in the original inverse 
seesaw formulation [23], the smallness of fix was attributed to the supersymmetry breaking effects 
in a superstring inspired Eq scenario. It is also possible to keep fix small by relating it to a tiny 
VEV generated dynamically. An analysis along this line was carried out in an extended version of 
the NMSSM [29J, where the origin of small VEV can be traced to the assumption of a vanishing 
trilinear scalar coupling at the GUT scale. A small fix can also be realized in a supersymmetric 
SO(10) context [30] . Regarding the second question, we recall the example of a non-supersymmetric 
SO(10) framework, which also contains the remnants of a larger Eg symmetry [27], where fix is 
generated at two-loop level, but /xjv is generated at a higher loop justifying its relative smallness. 
In our work, we do not advocate any specific GUT scenario to provide the dynamics that generates 
a small fix- We treat fix as an effective phenomenological parameter of unspecified origin, whose 
smallness derives its origin in some unknown hidden sector dynamics. We simply set its value to 
reproduce the correct light neutrino mass. 

We now illustrate the pattern of neutrino masses with only one generation. In the {^f v , ^n, ^x} 
basis, the (3 x 3) neutrino mass matrix is given by 

/Or, 

M = \ m D M N I , (6) 






m D 





m D 





M N 





M N 


fix 



yielding the mass eigenvalues (mi <C ttt-2, 
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m W _ _ ^ /, f 2 I _2 , M N^X 



mi = ml + Ml > m ^ = ^ M " + mD + 2{ml + Ml) - (7) 

The important thing to observe here is that the lightness of the smallest eigenvalue mi is due to 
the smallness of fix- The other two eigenvalues (m,2 or 7713) can have a mass around 10 GeV, and 
their presence significantly influences the decay pattern of A\. 

We now compute the branching ratios of Ai into the invisible modes comprising of the ty u , 
and the *&x states. Rigorously speaking, one should first diagonalize the mass matrix of Eq. ([6]) 
to determine the physical neutrino states. However, for our purpose it suffices to estimate the 
branching fractions of Ai into the ^ u ^n and ^jv^x interaction states. Recall from Eqs. ([I]) and 
([5]) that the decay of Ai into ^ U ^N will depend on how large the doublet component of Ai is, 
i.e. on how large cos 6 a is, whereas the decay into ^n^x will depend on the amount of As inside 
Ai, i.e. on the magnitude of sin^. Below, we present the branching ratios into invisible modes 
normalised to the visible ones (neglecting, for simplicity, the phase space effects). 



Br (Ai ->■ ^ y ^ N ) m 2 D 



Br (Ai ->//)+ Br (Ai -> cc) m 2 f tan 4 (3 + m 2 c ' 



(8) 



Br(Ai^^x) „ ^ 9a __^__* (9) 



Br (Ai ->//)+ Br {A x -> cc) m 2 f tan 2 + m 2 c cot 2 v 2 ' 
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where v = J v\ + v\ ~ 174 GeV. Notice that the dominant visible decay modes of A\ are //(/ = 
H,T,b) and cc. Of course, the cc mode would be numerically relevant if itia 1 < Ira^ and tan/3 is 
small. Note that the branching ratio into ^n^x dominates over that into for two reasons 

- firstly, there is a tan 2 6 a prefactor for the former which can be rather large if A\ has a dominant 
singlet component; secondly, if the m 2 term in the denominator of the branching ratio expressions 
is numerically relevant, then the 'Iv^'yv channel suffers a suppression by an additional tan 2 /3 factor. 

For a numerical illustration, we make two choices of tan/3 = (3,20), and fix cos 9a = 0.1, which 
yield = cos 9a tan (3 = (0.3,2). We recall that the upper limit on for rriA 1 < 8 GeV in the 
minimal NMSSM has been obtained primarily from radiative T-decays, and the limit is between 
0.7 to 3.0 for tan/3 = 50, while it is 30 or above for tan/3 = 1.5 [31]. A value of X^ = 2 is in 
fact slightly above the upper limit for mA 1 in the range of 4 to 8 GeV. In the present scenario, 
A± has a significant branching ratio into invisible modes which, in turn, considerably relax the 
upper bound on X&. Here we do not choose a very large value of tan /3 as that would increase the 
branching ratio of A\ into visible modes. The value of m,A 1 is chosen to be somewhat larger than 

1 /2 

Mat, so that the phase space suppression, given by the factor N 1 — (^p) 2 j / jl — (Ir/^) 2 ' 
is not numerically significant. We consider two values for = (5, 30) GeV. The rationale behind 
choosing M/v = 5 GeV is that it allows us to explore rriA 1 < 10 GeV, a regime where constraints 
from T- and -B-decays are particularly restrictive - see Table [TJ On the other hand, the choice 
Mjv = 30 GeV implies that A\ is moderately heavy (mA L > 30 GeV) which corresponds to the 
range where LEP and -B-decay constraints are relevant. We display our results in Table [2j For 
numerical illustration, we have assumed vs ~ 0(v). The main conclusion is that if cos 6a is small, 
A± has a dominant singlet component (which is generally the case when A\ is light [4]), then for a 
reasonable part of the parameter space A\ can have a sizable invisible branching ratio which would 
weaken many of the constraints discussed in the beginning. However, it is important to stress that 
cos 6a should not be excessively small, since in that case the purely singlet A\ would be completely 
decoupled from the visible sector. 





tan/3 = 20, cos(9 A = 0.1 


tan/3 = 3, cos 6 a = 0.1 


M N (GeV) 


5 


30 


5 


30 


Br (Ai -> 9 V 9 N ) 


7 x 10~ 5 


3 x 10~ e 


4 x lO" 3 


1 x 10" 4 


Br {A x -> * N *x) 


0.7 


0.9 


~ 1 


~ 1 



Table 2: Invisible branching ratios of the lightest NMSSM pseudoscalar for mo = 10 GeV, Mn = (5,30) 
GeV, and fix = 1 eV. 

What about the CP-even Higgs mass limit? Since the ZZh coupling is diluted with respect 
to its SM value as a result of singlet admixture, the direct search (lower) limit on from its 
non-observation will be lower than the SM limit of 114 GeV. A study based on the OPAL data 
from the LEP-2 run shows how the lower limit on decreases (assuming the Higgs production via 
Higgs-strahlung process) as £ = a(Zh)Bi(h — > invisible) /a SM (Zh) gets smaller than unity [32]. In 
our case, we have not only a mixing between the MSSM part of the CP-even Higgs and the singlet 
CP-even component, but also a sizable branching ratio of h — > A\A\ — > invisible. The lower limit 
on the lightest CP-even Higgs mass will then decrease accordingly. 

To conclude, we explore the possibility of having a very light (of order 10 GeV) pseudoscalar in the 
NMSSM, which has interesting consequences for CP-even neutral Higgs search at colliders, as well 
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as for facilitating dark matter annihilation. In this context, we have extended the minimal NMSSM 
with two additional gauge singlets carrying opposite lepton numbers for two specific reasons. On 
one hand, they provide a substantial invisible decay channel to the lightest pseudoscalar which 
helps relaxing or even evading some of the tight constraints from T- and U-decays. On the other 
hand, they naturally set up the stage for implementing the inverse seesaw mechanism in order to 
generate light neutrino masses. What is phenomenologically interesting is that this can be done 
using order one neutrino Yukawa couplings and employing neutrino Dirac masses of a few tens of 
GeV. To account for the experimental values of the two mass squared differences and the three 
mixing angles of light neutrinos, one would of course have to extend the number of N and X 
superfields. 
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